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ABSTRACT 

Modelling of high resolution Balmer line profiles in the early-time spectra of 
SN 1998S shows that the inferred fast (w 400 km s~^) circumstellar (CS) gas on 
days 23 and 42 post-explosion is confined to a narrow, negative velocity gradient shell 
just above the photosphere. This gas may be identified with a slow {v < 40 km s^^) 
progenitor wind accelerated at the ejecta-wind interface. In this scenario, the photo- 
sphere coincides with a cool dense shell formed in the reverse shock. Acceleration by 
radiation from the supernova or by a shock-accelerated relativistic particle precursor 
are both possible explanations for the observed fast CS gas. An alternative, equally 
plausible scenario is that the fast CS gas is accelerated within shocked clouds engulfed 
by the outer shock, as it propagates through the intercloud wind. 

Key words: supernovae - circumstellar matter - : stars. 



1 INTRODUCTION 

Supernova (SN) 1998S in NGC 3877 was discovered by Li & 
Wan (1998), probably before maximum light (Fassia et al. 

2000) . The presence of narrow lines in its early-time spec- 
tra led to its classification as a type Iln (SN Iln) (Leonard 
et al. 2000). Nevertheless, the close resemblance of its light 
curve and spectral development to those of SN 1979C (Liu 
et al. 2000) means that it is also valid to classify SN 1998S 
as another bright SN IIL. Whatever the assigned class, there 
is strong evidence for the existence of a dense circumstellar 
(CS) wind in its immediate vicinity of SN 1998S. This in- 
cludes the presence of narrow CS lines (Bowen et al. 2000; 
Fassia et al. 2001), broad emission lines originating from the 
interaction (Leonard et al. 2000; Gerardy et al. 2000; Fas- 
sia et al. 2001), and radio and X-ray emission (Pooley et al. 

2001) . Despite extensive effort in the past, we still have a 
poor understanding of the origin of the large pre-supernova 
mass loss and the ejecta-wind interaction phenomenon in 
SNe Tin. 

SN 1998S is the most comprehensively observed SN Iln, 
thus providing us with an excellent opportunity for gaining 
a deeper insight into the morphology and origin of the dense 



CS environment and the physics of the interaction. Of par- 
ticular interest is the origin of two distinct CS wind velocity 
components deduced from high resolution early-time spec- 
tra of SN 1998S taken at ~2-5 weeks after the explosion 
(Bowen et al. 2000; Fassia et al. 2001). These had velocities 
of Uw = 40 — 50 km s~^, and u ~ 400 km s~^, respectively. 
The slower component is revealed by both narrow H I and 
He I P Cygni profiles and forbidden lines, while the faster 
component is seen only in broader H I and He I absorp- 
tion lines. The minima of the broad Ha and H/3 absorptions 
shifted from ~-115 km s"^ on 1998 March 20 to ~-85 km 
s~^ by 1998 April 8. Similarly broad CS absorption lines 
are also seen in high resolution ultraviolet (UV) spectra via 
the Mg II 2800 A doublet (Bowen et al. 2000). While the 
presence of the slow wind is unsurprising for a SN II progen- 
itor, the origin of the faster CS material (CSM) is not clear. 
Bowen et al. (2000) suggest it originated as a pre-existing 
fast wind component (fast shell). Fassia et al. (2001) suggest 
it was due to a faster wind produced during a recent BSC 
phase. A third possibility is that it was due to pre-shock 
radiative acceleration of the slow wind (Lentz et al. 2001; 
Fassia et al. 2001). 

In this paper we shall focus on the nature of the fast 
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CS component which, we beheve, is related to the structure 
of both the wind and the ejecta-wind interface. We start 
with a study of the kinematics and location of the fast CS 
gas (Section 2), making use of the high resolution profiles of 
Balmer lines observed on 1998 March 20 and 1998 April 8 
(Bowen et al. 2000; Fassia et al. 2001). Through our mod- 
elling of the CS gas line profiles, we reveal that the fast CS 
component appears to lie very close to the SN photosphere. 
Indeed, this proximity presents a challenge in understanding 
the structure of the ejecta-wind interface. A solution will be 
proposed. We then consider the plausibility that the fast CS 
gas results from the pre-shock acceleration of the slow wind 
by radiation or relativistic particles from the SN. (Section 3). 
As an alternative explanation, we also investigate shocked 
dense wind clumps (Section 4). We note that the CS wind 
considered here corresponds to the 'Mid-CSM/Outer-CSM' 
regions defined in Fassia et al. (2001). This will be discussed 
further in Section 5, where we conclude that the division into 
Mid- and Outer-CSM zones is probably not required after 
all. The 'Inner-CSM' region (Fassia et al. 2001), overrun by 
the ejecta around 2 weeks post-explosion, is examined by 
Chugai (2001). 

Based on Chugai's (2001) model for the SN-CSM inter- 
action and the light curve, we prefer 1998 February 25 as the 
date of the SN 1998S explosion. Adopting this as zero epoch, 
March 20 and April 8 become days 23 and 42, respectively. 
This should be kept in mind when comparing this work with 
other papers on SN 1998S. 



2 MODELLING THE CIRCUMSTELLAR LINE 
PROFILES 

In the standard picture of spherically symmetric SN ejecta 
interaction with a smooth CS wind, the ejecta-wind interface 
has a two-shock structure with the outer and inner shocks 
running outward into the wind and inward into the ejecta, 
respectively (Chevalier 1982; Nadyozhin 1981, 1985). We de- 
fine a wind density parameter (linear density) w = 47rr^p, 
(g cm~^), where p is the density of the unshocked CS wind. 
With w ~ 10^^ g cm~^ , which is typical for SNe II with dense 
winds (e.g. SN 1979C, see Lundqvist & Fransson 1988), the 
inner shock wave is radiative. Consequently, a cool (T ~ 10* 
K) dense shell (CDS) of shocked SN ejecta forms between the 
outer and inner shocks (Chevalier & Fransson 1985, 1994). 

To investigate the origin of the CS line profiles, we con- 
sider a model of the line-forming region specified in terms 
of a) the photospheric radius, Rp, b) the extent of a shell of 
fast CS gas, Ri < r < R2, and c) the extent of a slow exter- 
nal wind, R2 < r < R3. We express the radii in terms of the 
photospheric radius, i.e. Rp = 1, unless the physical units 
are needed. The model is similar, in some respects, to that 
of Lentz et al. (2001), but with two important diff'erences. 
First, in their model there is a significant "SN atmosphere" 
zone lying between the SN photosphere and fast CS gas 
(-Rp < r < _Ri). In our model the SN atmosphere is absent 
(we explain why below). There may be a layer of very hot 
ionised post-shock gas between the photosphere and the fast 
CS gas {Rp < r < Ri), which does not contribute to line 
emission or absorption in the optical, but reprocesses the 
photospheric spectrum in the ultraviolet due to Comptoniza- 
tion (Lundqvist et al., in preparation). Second, in Lentz et 
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Figure 1. High resolution Ho spectra of SN 1998S for two epochs 
(Fassia et al. 2001). The flux is in units of 10~^* erg cm~^ 
A~^. The lower spectrum is shifted down by the value shown in 
parentheses. The symbols a-e show features referred to in the 
text. 



al. the velocity within the fast CS gas is assumed constant, 
while in our model the velocity can vary with radius. Since 
the kinematic structure of the resulting fast CS gas is not 
known, we consider initially three different possibilities: (a) 
a fast, constant velocity wind, (b) an accelerated fiow with a 
positive velocity gradient, i; oc r, (this might arise from a CS 
shell lost in a violent mass ejection before the SN explosion 
(Grasberg & Nadyozhin 1987)), (c) an accelerated flow with 
a negative velocity gradient {dv/dr < 0), characteristic of 
pre-shock acceleration of the slow wind. For the case with 
a negative velocity gradient, we adopt as the description of 
the fast CS gas velocity u: 

V -J 



= {U[ — Uw) 



■.R2-R1 



(1) 



where u„, the velocity of the slow wind, is assumed to be 
constant with radius. Thus the velocity of the accelerated 
CS gas is Uf at its innermost limit _Ri, falling to «w at its 
outermost limit R2 where it merges with the unaccelerated 
CS wind. Note that u; <^ 1 corresponds to a very steep 
velocity gradient, while ui ^ 1 gives a mild slope. 

Given the population of the second level n2{r), and the 
line source function S{r), the CS gas line profile can then 
be easily calculated for the velocity fields specified in (a)- 
(c) above. We do not compute the actual hydrogen level 
population in the wind (although this is itself an interesting 
problem). Instead, we parameterize the behaviour of n2(r) 
and S{r). Specifically, we adopt a power law distribution for 
the second level population n2(r) oc for r > R\, with 
712 (-Ri) and p as free parameters. We split the Ha source 
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Table 1. Parameters of illustrative fast CS gas line profile models 
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function into net emission and scattering terms, respectively 
Se and Sac, with Sc oc r"' and Ssc = W'-^'c, where W is the 
dilution factor and Jc is the continuum intensity. This repre- 
sentation reduces uncertainty in the Ha source function. For 
other Balmer lines we cannot assume such a simple represen- 
tation of the source function. However, the scattering term 
contribution is expected to be substantially lower than for 
Hq, and so we simply adopt the power law S oc r~' for these 
lines. Rather than using the Sobolev escape approximation, 
we directly integrated the radiative transfer equation, adopt- 
ing a Gaussian local profile with a thermal most probable 
velocity Wth = ^IkT /m. The model parameters comprised 
R\, R2, R3, Wf, Ww, Mth, a;, p, q, n2{R\) and S{Ri). 

To assess the viability of the different CS kinematics 
suggested above we shall concentrate on the kinematics of 
the fast flow, but always including the slow wind compo- 
nent in the model. In all the models the maximum velocity 
of the fast flow, U{ = 500 km s~^. The kinematics of the 
slow wind is assumed to be a combination of its steady out- 
flow velocity ii„ — 30 km s^^ within which the hydrogen 
gas has a thermal most probable velocity uth — f5 km s~^ 
(corresponding to ~ 14000 K). The adopted values for 
and Wth are based on our best fits to the observations for the 
slow wind lines (see below), although their precise values are 
not crucial here. A successful model of the fast CS line com- 
ponent needs to be able to reproduce the specific features 
indicated in Fig. f , viz. (i) the low velocity of the absorption 
minimum (point b); (ii) the slowing down of the absorbing 
gas velocity between day 23 and day 42, as indicated by 
the blue wing (a-b); (iii) the relative weakness of the red 
emission wing (e-f), presumably indicating an occultation 
effect. 

The resulting model profiles are displayed in Fig. 2 with 
the parameters given in Table f . Comparison of Figures f 
and 2 shows that model f (fast, constant velocity wind) pre- 
dicts a velocity for the absorption minimum which is far 
greater than that observed. Variation of the free parame- 
ters (J?i, R2, p, q) does not improve the fit. We therefore 
rule out this model. A somewhat better match is achieved 
with model 2 {v (x r). Moreover, if we increase the outer 
radius by factor of three, then since we have a fixed max- 
imal velocity, the velocity of the absorption minimum can 
be made even lower. However, this model predicts that the 
characteristic velocity derived from the fast CS gas line pro- 
file should increase with time as SN overtakes the innermost, 
lower velocity CS gas. In fact we observe that the charac- 
teristic velocity decreases with time (Fig. f), and so we also 
discard this model. 

Thus, we are left with the negative velocity gradient 
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Figure 2. Line profiles predicted by a two-component model 
consisting of an inner fast flow and an outer slow wind. The upper 
panel shows the constant velocity fast wind Model f (solid line) 
and the v oc r Model 2 (dotted line). The middle and lower panels 
show negative velocity gradient model profiles. The middle panel 
shows the best-fit Model 3 (thick line), compared with Model 4 
which has a larger outer radius R2 = f.5 (thinnest line), and 
Model 5 which has a gap between the photosphere and the inner 
radius (Ri = 1.1). The lower panel also shows Model 3 (thick 
line), compared here with Models 6 (thinnest line) and 7 which 
have, respectively, steeper (oj = f.3) and milder (oj = 0.3) velocity 
gradients. 

class of models (Models 3-7). In Fig. 2 (middle panel) we 
show Model 3 which gives a good fit to the observed day 23 
profile (see below). With Models 4-7, we illustrate the sensi- 
tivity of these negative gradient models to parameter varia- 
tions. Increasing the outer limit, R2 of the accelerated CSM 
(Model 4) (Fig. 2 middle panel) produces a red wing which 
is much stronger than that observed. The introduction of a 
gap, SR = Ri — Rp = O.I, between the photosphere and the 
fast CS gas (Model 5) also results in an unacceptably strong 
red wing (Fig. 2 middle panel). Moreover, both Models 4 
& 5 exhibit significant differences from Model 3 at radial 
velocities of ~ —100 km s^^. Variation of the velocity gradi- 
ent index, cu, also has a strong effect. This is illustrated by 
models 6 & 7, again compared with Model 3 (Fig. 2, lower 
panel). Thus, given a judicious selection of parameters, these 
preliminary calculations demonstrate that a model compris- 
ing a fast CS shell with a negative velocity gradient can 
successfully reproduce the observations. 

We now fit the Ha CS line profiles for both epochs and 
compute the other Balmer line profiles using the same pa- 
rameters as for Ha, tuning only the amplitude of the source 
function. A redshift of -1-850 km s~^ (Fassia et al. 2001) 
was adopted. We explored the model parameter space and 
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Figure 3. High resolution Ha (top panel) and H/3 spectral pro- 
files in SN 1998S on day 23. Thin line - observation, thick line - 
model (Table 2) 



found that the observed Ha profile could only be reproduced 
with a narrow (or zero) gap between the photosphere and 
the fast CS gas, viz. Ri < 1.057?p, with Ri — Rp being 
the best choice. The power-law indices, p, q and amplitudes, 
^2(^1), S{Ri) of, respectively, the second level population 
and source function were estimated by imposing the require- 
ment that both the fast and slow CS absorption components 
of all the hydrogen lines should be reproduced. The slow 
wind outflow velocity Mw and thermal velocity itth were fixed 
by the positions and widths of the narrow lines. In general, 
the slow wind outer radius, R3, is poorly determined. We 
tested values in the range 2.5 < 7?3 < 8 but were always 
able to produce acceptable fits by adjusting p and q. We 
have, therefore, arbitrarily adopted .Ra = 4 on day 23 and 
J?3 = 3 on day 42. We judge that the Balmer line profile 
fits for days 23 and 42 (Figs. 3 & 4) are quite good. The 
parameters used are given in Table 2. In Fassia et al. (2001) 
we found that the narrow H-line widths corresponded to a 
velocity of ~60 km/s, and we suggested that part of this 
was due thermal broadening. This is quantified here by our 
modelling of the narrow hydrogen lines which indicates a 
slow wind velocity of 30 km s~^ for hydrogen. Curiously, 
this is actually slightly slower than the ~40 km/s wind ve- 
locity inferred from the narrow forbidden lines (Fassia et al. 
2001). 

The model of the fast CS gas defined by the parame- 
ters in Table 2 is henceforth dubbed the 'kinematic model'. 
The essential feature of this model is that the fast CS gas is 



confined to a narrow layer, AR/R ~ 0.3, nearly attached to 
the photosphere. This indicates that the photosphere resides 
in the outermost layers of the SN ejecta and possibly coin- 
cides with the CDS. The latter point, in fact, agrees with 
recent modelling of the SN 1998S interaction with the CS 
gas, which demonstrates that the CDS remains opaque in 
the Paschen continuum for about 45 days (Chugai 2001). 
The essential elements of the kinematic model of the fast 
CS gas are illustrated in Fig. 5. 

This picture is, however, somewhat disturbing since it 
leaves no place for the outer shock wave. The point is that 
the fast CS gas presumably lies in the pre-shock zone and 
therefore must be separated from the CDS by the post-shock 
layer. The extent of this layer depends on whether the shock 
is adiabatic or radiative. In the adiabatic case the typical 
thickness of the postshock layer for the outer shock is about 
AR/R ^ 0.2 (Chevalier 1982). This is in conflict with the 
kinematic model of the pre-shock fast CS gas, which re- 
quires AR/R < 0.05. On the other hand, if the outer shock 
were essentially radiative, the post-shock layer might be thin 
enough to meet this requirement. However, the problem with 
the latter possibility is that the wind density in the region 
r > 10^^ cm, where the interaction proceeds after March 16 
(day 19) is very low, w ~ 10^^ g cm~^, according to the anal- 
ysis of the fast CS ultraviolet lines on March 16 (Lentz et al. 
2001). (We note also, that this is consistent with the lower 
limit of ~ 10^* g cm~^ on day 42 which can be obtained 
from the [O III]-derived density lower limit and the shock 
radius (Fassia et al. 2000, 2001)). For such a wind the outer 
shock is essentially adiabatic. We therefore face a 'proximity 
problem', namely, how to reconcile the proximity of the fast 
CS gas to the photosphere on days 23 and 42 with the claim 
that the outer shock wave is adiabatic at these epochs. 

Two solutions to this problem are conceivable. The first 
possibility is that the outer shock is adiabatic, i.e R\ m 1.2, 
but that the CDS is strongly corrugated, presumably due 
to Ray leigh- Taylor instabilities, in such a way that spikes 
of dense opaque CDS material protrude into the postshock 
layer which lies in the region _Rp < r < 7?i . The morphology 
of the corrugated CDS in this case is similar to the 'hedge- 
hog' model used by Cid Fernandes & Terlevich (1994) to 
explain line profiles in SNe II interacting with a dense CS 
environment. The thick forest of spikes can, in principle, 
produce an efficient limb occultation of the far side of the 
fast CS shell and thus bring about the effect of fast CS gas- 
photosphere proximity. 

The second solution is based upon the concept of a 
clumpy wind. Wind dumpiness is a rather common phe- 
nomenon for AGB and post-AGB stars (Olofsson 1994) and 
is occasionally revealed in interacting SNe II (Chugai 1993). 
We envisage a two-phase wind with dense clouds embedded 
in the rarefied intercloud gas. The SN interaction with this 
wind results in a fast outer shock running through the in- 
tercloud gas, and slow shocks running into the dense clouds. 
Thus, the nearly-stationary clouds become engulfed by the 
outer shock and, for a certain density contrast and cloud 
size, can survive in the postshock gas, eventually encounter- 
ing the CDS. If we attribute the fast CS gas to the effect 
of shocks running into the engulfed clouds then the proxim- 
ity problem may be resolved. Two options are conceivable, 
depending on whether the cloud shock is adiabatic or radia- 
tive. In the first option, the mechanism which accelerates 
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Figure 4. High resolution Balmer line profiles from Ho (top 
panel) to H5 (bottom panel) in SN 1998S on day 42. Thin line - 
observation, thick line - model (Table 1). 



Table 2. Parameters of 'best' fit models 
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the unshocked slow gas in the clouds up to the fast CS gas 
velocity is presumably the same as would be the case for a 
smooth wind in the pre-shock region. In the case of radia- 
tive cloud shocks the fast CS gas may be identified directly 
with the shocked cloud gas. The latter possibility will be 
considered in detail in Section 4. 



3 FAST CS GAS ARISING FROM PRE-SHOCK 
ACCELERATION 

Between the two epochs (23 and 42 d) the photospheric ra- 
dius expanded by a factor of ~ 1.5 (Fassia et al. 2000). Yet 
the above analysis indicates that, at both epochs, the fast 
CS gas resided in a thin shell (AR/R ~ 0.3) attached to the 
photosphere. Moreover, the gas flow in this shell has a neg- 



ative velocity gradient, declining from the maximal velocity 
to ~ 30 km s~^ at the inner edge of the disturbed slow wind. 
All this leads us to the proposition that the fast CS gas is 
not a pre-existing flow, nor is it accelerated shortly after 
the explosion by the radiation from the shock breakout, but 
instead is created by the ongoing acceleration of the slow 
wind in the immediate vicinity of the photosphere. The ac- 
celeration might be produced by the SN radiation (Chevalier 
1981) and/or a relativistic particle precursor (Chugai 1997). 



3.1 Radiative acceleration 

The acceleration of an ionised, optically thin wind via Thom- 
son scattering was invoked to explain the fast CS features 
observed in SN 1993J at early epochs (Fransson, Lundqvist 
& Chevalier 1996) . Ignoring the initial wind velocity, Thom- 
son scattering can accelerate the CS gas at radius r up to 
V = 80E4gr'[^ km s~^, where £49 is the time-integrated ra- 
diation energy in units of 10*® erg, and ris is r in units of 
10^^ cm. Fransson et al. (1996) estimated that £^49 ~ 1 for 
SN 1993J using explosion models and calculating the emis- 
sion from the circumstellar shock. For SN 1998S we have no 
information on the total luminosity for the early evolution 
until the estimates for the bolometric luminosity of Fassia 
et al. (2000), starting on day 17. Integrating the estimated 
bolometric luminosity between days 17 and 42, we obtain 
£49 ~ 7 if we assume that the spectrum is blackbody. There 
is probably considerable uncertainty in this estimate. We 
have not included the early emission from the photosphere 
and from the circumstellar interaction, which may lead to 
an underestimate, but we have also ignored the effects of 
UV line blanketing in the blackbody approximation, thus 
producing a possible overestimate. However, the bolomet- 
ric luminosity is certainly considerably higher than for SN 
1993J, with a likely range of i?49 = 5—10. Remarkably, in 
the SN II models of Falk & Arnett (1977) with an extended 
presupernova envelope of a size and density similar to those 
of the SN 1998S inner sheU (Chugai 2001), the radiated en- 
ergy S49 ~ 10. 

On day 42, ns ~ 2 (Fassia et al. 2000) which, with 
Eia = 5- 10, yields « ~ 100 - 200 km s~^ for the CS gas. 
This is faster than the slow wind velocity of ~ 30 km s~^, 
suggesting that the observed slow wind may have undergone 
some radiative acceleration following the supernova explo- 
sion. Yet these velocities are insufficient to account for the 
fast CS shell at this epoch. The same is true for day 23 when 
ri5 ~ 1.3. The value of £^49 is even more uncertain at that 
epoch, probably in the range E^g = 1 — 6, using the results 
of Fassia et al. (2000) . The corresponding velocity range for 
the CS gas is then 1; ~ 50 — 300 km s~^, again too low for 
the fast CS gas. We therefore have to seek an additional ac- 
celeration mechanism to account for the velocity of the fast 
CS gas. 

The radiative force can be enhanced by including the 
effects of UV line absorption. This mechanism was exploited 
by Castor, Abbott & Klein (1975) to account for hot star 
winds. It was also invoked by Chevalier (1981) to explain the 
fast CS features in SN 1979C via radiative acceleration in 
the pre-shock zone, and discussed by Fransson (1986) who 
argued that lines could pre-accelerate the gas up to x4 — 10 
more efficiently than Thomson scattering. 
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Figure 5. Illustration of the principal components of the kine- 
matic model for the fast and slow CS line formation. 



The radiative force per unit volume due to the line ab- 
sorption effect (Castor et al. 1975) is 



f 



dr 



(2) 



where Fi, is the monochromatic flux. Tj is the radial Sobolev 
optical depth in the line, where j is summed over the total 
number of lines with the lines arranged in order of frequency. 
From the dependence of the velocity gradient on / we see 
that line acceleration is a runaway effect; the gas that has 
been accelerated is also the gas that most easily gets fur- 
ther accelerated. This means that a steeper velocity profile 
can be created than with pure Thomson scattering (see also 
Fransson 1986), thus providing an explanation for the line 
profiles in Section 2. Another mechanism which can steepen 
the velocity profile is related to the degradation of the UV 
radiation due to photon splitting. This results in a more ef- 
ficient acceleration in the inner layers of the pre-shock zone, 
where the UV fiux is greater. Clearly, an adequate treat- 
ment of the line effects in radiative acceleration requires the 
solution of the radiative transfer with the statistical equi- 
librium of ions and level populations. This is a formidably 
complicated task, beyond the scope of this work. We shall, 
instead, restrict ourselves to estimating the radiative force 
and then the velocity increase due to the lines effect. 

Let ^ be the ratio of line absorption force / to the Thom- 
son scattering force /t = aTn^F/c {F is the total radiation 
fiux integrated over wavelength). Using the equation of mo- 
tion, the velocity change of the fast CS gas due to the lines 
effect in the pre-shock zone AR may then be written ap- 
proximately as 



(3) 



where R is the radius of the outer shock, At = 0.3 cm^ 
g~^ is the Thomson opacity, and Vb is the outer shock wave 
velocity, which we also identify with the velocity of the CDS. 
To estimate the radiative force due to line absorption 



we use the list of 1.5 x 10^ lines with log [gf) > —5 com- 
piled by Ron Eastman (Blinnikov et al. 1998) and extended 
here using the R.L. Kurucz Data Base available on the 
web. Solar composition and a Saha-Boltzmann distribution 
are assumed, while the photospheric radiation is treated 
as black-body with the flux = TiB^iT). In Fig. 6 we 
show C, = ///t as a function of temperature for densi- 
ties of 2 X 10"", 2 X 10"" and 2 x 10"^^ g cm"^ (We 
note that these values are consistent with the lower limit 
Tie > 10"^** cm~^ obtained by Fassia et al. (2001) from 
observed [O III] line ratios.) The velocity gradient in Eq. 
(|) is taken to be Av/AR, with Av = 400 km s"-^ and 
AR = 3 X lO'^* cm. The results indicate that is lower for 
higher densities. This is due to line saturation effects [see 
Eq. (^)]. The non-monotonic temperature behaviour is due 
to the strong temperature dependence of the ionisation frac- 
tions in different ions. The low temperature peak at 6000- 
7000 K is caused primarily by the low ionisation ions of Fe. 

On day 23, with the photospheric radius at 1.3 x 10^^ 
cm (Fassia et al. 2000), the average radius of the fast CS 
shell is « 1.4 x 10^^ cm. Therefore, a wind linear density 
w « 10^^ g cm"^ would correspond to the density 5 x 10"^^ 
g cm"^. With the equilibrium temperature ~ 1.4 x lO** K 
(Fassia et al. 2000) we have thus ^ « 1 (Fig. 6). Adopting 
Eq.(|) i? = 1.4 X lO^'^ cm, = 6000 km s"^ (Chugai 2001), 
luminosity L = 3 x 10'*'' erg s"^ (Fassia et al. 2000), and 
AR/R = 0.3, we obtain from Eq.(||) Au{ ^ 90 km s"\ The 
integrated acceleration due to Thomson scattering was esti- 
mated above as 50 — 300 km s"^ for this particular epoch. 
In combination with the lines effect, this yields 140 — 390 
km s"* I.e. the upper limit is about equal to the observed 
maximal velocity of the fast CS gas (~ 400 km s"^). Thus, 
radiative acceleration is a plausible mechanism for the ac- 
celeration of the fast CS gas. This possibility is enhanced 
when we note that the effect of the lines may have been 
underestimated since (a) the line list may be incomplete, 
(b) multiple scattering of radiation in the UV line forest 
has been ignored, and (c) the photospheric spectrum is not 
strictly black-body. Also, as has been shown by Lundqvist 
& Fransson (1988) and Fransson et al. (1996), the wind con- 
ditions can be far from thermodynamic equilibrium. Most 
likely it is 'over-ionised' which would make the wind behave 
like a gas which was hotter (c/. Fig. 6) than the 1.4 x 10'* K 
we have assumed, thereby increasing the radiation force. 

The slow-down of the fast CS gas between days 23 and 
42 can be naturally explained using Eq.(^. Using informa- 
tion from Fassia et al. (2000), the ratio L/Rvs decreases by 
a factor > 2 during this period, and unless increases signif- 
icantly during the same period, the line acceleration shows 
roughly the right temporal behaviour. 



3.2 Relativistic particles 

An alternative mechanism for producing pre-shock accelera- 
tion is the dynamical effect of the expected relativistic par- 
ticle shock precursor (Kang & Jones 1997). Let a fraction rj 
of the kinetic luminosity, {1/2)wVb, dissipated in the outer 
shock, be converted into relativistic particles trapped in the 
shock wave region. In the approximation of constant ejecta 
velocity, Vs, the energy equation with adiabatic losses gives 
as the total energy of relativistic particles Er = {1/ A)r)RwvB 
(where R is the shock wave radius) , and the relativistic par- 
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Figure 6. Ratio of the cumulative line absorption radiation force 
to the Thomson scattering radiation force, as a function of the 
temperature. Three lines correspond to densities p = 2 X lO"^'^ g 
cm~'^ (dashed), p = 2 X 10~^^ g cm~^ (solid), and p = 2 X 10^^^ 
g cm~'^ (dotted). 

tide pressure in the precursor Ap = (l/3)^-E'r/(47rJ?^A_R) 
(where ^ is the fraction of relativistic particles in the pre- 
cursor). The equation of motion of the pre-shock flow with 
velocity «i in the shock frame, is pvi = —Ap/AR, and so 
we have as the velocity of the pre-shock gas in the lab frame 

""^^h (^) ^''^^ 

Setting AR/R = 0.3 and «s = 6000 km s~^, to obtain a 
precursor velocity of ut = 500 km s~^ we need ^t] ~ 0.3. 
Given comparable fractions of relativistic particle energy in 
the post-shock and pre-shock regions (Kang & Jones 1997), 
we may adopt ^ ~ 0.5. Thus, to accelerate pre-shock matter 
up to 500 km s~^, the conversion efficiency of shock energy 
into relativistic particles must be ~ 0.6. While this value is 
close to the absolute upper limit, it is nevertheless plausible 
given the high efficiency of particle acceleration in high Mach 
number astrophysical shocks (Blandford 1984). Thus, from 
energy considerations alone, we cannot rule out the possibil- 
ity that the relativistic particle precursor may significantly 
contribute to the acceleration of the wind in SN 1998S. 



4 FAST CS GAS ARISING FROM SHOCKED 
WIND CLUMPS 

We now consider a radically different mechanism for pro- 
ducing the fast CS gas. This has already been mentioned in 
Section 2. It essentially exploits the possibility that the wind 
may have a clumpy structure. In this scenario, the fast CS 
gas is identified with radiatively-cooled shocked gas in wind 
clouds engulfed by the intercloud shock. This may, in princi- 
ple, reproduce both the characteristic velocity of the fast CS 
gas and resolve the photosphere proximity problem. More- 
over, it may also produce a negative velocity gradient. Since 
the postshock density Pps increases from the outer shock 
jump towards the contact discontinuity (Chevalier 1982), it 



means that the cloud shock velocity Vc ~ i's(pps/ Pc)~^^^ also 
increases inwardly. Thus, the negative velocity gradient may 
well be mimicked by the ensemble of cloud shocks. 

A specific feature of the proposed mechanism is a 'strip- 
ping flow' produced by the mass loss fi-om the shocked 
cloud due to a Kelvin-Helmholtz instability (Klein, McKee 
& Colella 1991). The stripping flow is characterized by a 
broad spectrum of fragment sizes and velocities with the 
smaller fragments acquiring the higher velocities. Tiny, high 
velocity fragments are subjected to efficient mixing with the 
hot outer shock gas, which results in their heating and ioni- 
sation. Therefore, we expect that the high velocity stripping 
flow should be revealed particularly by the presence of high 
ionisation lines. In fact HST UV spectra do show high ve- 
locity (several xlO^ km s~^) blue wings in C IV and Si IV 
resonance lines (Lundqvist et al., in preparation), and this 
may well be the consequence of such a stripping fiow. 

We now test the clumpy wind model against particu- 
lar observational requirements, particularly on day 23. The 
cloud/intercloud density ratio x must be consistent with the 
cloud shock wave velocity, «c = 2vsX~^^^ (the factor of two 
is from the factor of four density jump at the intercloud 
shock). With Vc ~ 500 km s~^ and a CDS velocity at this 
epoch of about Vb ~ 6000 km s~^ we thus obtain x ~ 600. 
To meet the requirement of proximity to the photosphere, 
the lifetime, tc = 2a/vc, of a cloud of radius a must be equal 
to or greater than the time, AR/v, which it takes to cross 
the intercloud post-shock layer [AR is the thickness of this 
layer). This leads to the constraint a/AR > 1/y/x ~ 0.04. 
We consider the case where the ratio is close to this lower 
limit, VIZ. a/AR ~ 0.05, in order to maximize the number 
of shocked clouds and so be consistent with the absence of 
strong fiuctuations in the line profile. Combining this ratio 
with AR/R ~ 0.2 and the value of the photospheric ra- 
dius at this epoch, R ^ 1.3 x W^^ cm (Fassia et al. 2000), 
we get a cloud radius a ~ 10^^ cm. The filling factor of 
the clouds (/) may be estimated from the condition that 
the covering factor of the shocked clouds is about unity, 
i.e. (3/4)/(A_R/a) ^ 1. From this, we obtain / ^ 0.07. 
Adopting all the cloud parameters derived above, we find 
that the number of clouds in the outer postshock layer is 
= ?.f(R/af(AR/R) 9 x lO''. 

We now demonstrate that, given the above choice of 
cloud parameters (a, Xi ^-nd /), the cloud shocks are es- 
sentially radiative. The condition that the cooling time is 
less than the time for the cloud shock wave propagation 
is a/uc > mp«c /(32Anc), where mp is the proton mass, 
and A ~ 2 X 10~^^ erg s~^ cm^ is the cooling function. 
This leads to the requirement one > 3 x 10^° cm~^, where 
Vc = 500 km s~^ is assumed. Given an average wind density 
n = ui/iinR^nip) « 2 x 10* cm"^, x = 600, and / = 0.07, 
we obtain a cloud density of ric ~ 3 x 10^ cm^^. Therefore, 
from the above lower limit for the cloud column density, aric, 
we find that for the cloud shocks to be radiative, a > 2 x 10^^ 
cm. This is consistent with the value a ~ 10^^ cm derived 
above using the observed values for Vc, Vs and R at 23 days. 

Finally, the adopted set of cloud parameters must be 
consistent with the day 23 observational upper limit for the 
net Ha luminosity L(Hq) < 10^^ erg s"^ (Fassia et al. 2001). 
To check this, we use the simple expression for the cumula- 
tive Ha luminosity of the shocked clouds 
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where ^li is the efficiency of Ha emission, of the order of 
10"^ (Chevalier & Fransson 1994). Substituting in Eq. (||) 
w = 10^*^ g cm-\ = 500 km s''^ , ip = 0.01, a/AR = 
0.05, one obtains L(Ha) < 10^* erg s^^, in accord with the 
observational upper limit. 

Thus, for the cloud shock scenario an appropriate set of 
cloud parameters might be x ~ 600, / ^ 0.1, a ^ 10^^ cm, 
with the average linear wind density w « 10^® g cm~^. Note, 
this model can account for the observed slow down with time 
of the average cloud shock velocity {vc ~ Va/^/x), if the 
density contrast x increases with radius between days 23 
and 42 and/or t^s decreases with time, both being plausible. 

Summing up, the above order of magnitude analysis 
shows that the shocked wind cloud model for the fast CS 
component may be matched to the bulk of the empirical 
constraints for a reasonable set of parameters. While the 
proposed model is quite crude, it does appear to offer a use- 
ful alternative explanation for the fast CS gas phenomenon, 
including the high velocity blue wings of the fast CS absorp- 
tion profiles in high ionisation UV lines. 



5 DISCUSSION 

The principal finding of this paper is that the fast CS gas ob- 
served in Balmer lines in SN 1998S during the era day 23-42 
has a negative velocity gradient flow confined to a narrow 
layer above, and attached to, the photosphere. The same 
should be true for the fast CS gas observed in He I (Bowen 
et al. 2000, Fassia et al 2001) and Mg II 2800 A (Bowen 
et al. 2000). We argue that the fast CS gas is the result of 
the acceleration of the slow wind in the immediate vicinity 
of the SN photosphere, and show that acceleration driven 
by radiation is a plausible mechanism. Thomson scattering 
alone, while signiflcant, is insufficient to account for the ob- 
served fast CS gas velocity. However, when we include the 
additional effects of UV line absorption, we can obtain the 
observed velocity. Moreover, this is a conservative conclu- 
sion given the that the line list may be incomplete, that we 
ignore multiple scattering in UV line forest, and that time- 
dependent and non-thermodynamic equilibrium effects are 
not included. We also demonstrate that the a relativistic 
particle precursor is able to accelerate pre-shock matter up 
to the observed fast CS gas velocity. 

The fast CS gas can also originate in a completely dif- 
ferent way; it could be the result of radiative shocks prop- 
agating into dense wind clumps together with the effects 
of a stripping flow. The stripping flow aspect of the clumpy 
wind model predicts the presence of very fast, highly ionised 
CS gas. This seems to be supported by the observation of 
high velocity blue absorption wings (several xlO'' km s^^) 
in C IV and Si IV (Lundqvist et al., in preparation). 

An intriguing fact, is that the narrow line width of the 
unaccelerated wind (~ 15 km s~^) indicates a rather moder- 
ate wind temperature ~ 10"* K. This is quite unlike the hot 
wind (~ 10^ K) at early epochs (t ~ 1 month) anticipated 
by the model of Lundqvist & Fransson (1988) for SN 1979C 
with comparable wind density, as well as for the somewhat 
less dense wind around SN 1993J (Fransson et al. 1996). 
This disagreement may indicate that either the outer shock 



is essentially radiative so that the Comptonization responsi- 
ble for the bulk of the UV radiation is inefficient, or that the 
wind is clumpy with a resulting low equilibrium temperature 
in the pre-shock wind clumps. We favour the second option, 
because a radiative outer shock would produce even higher 
UV luminosity than that of Comptonized radiation from an 
adiabatic shock. We note that a clumpy-wind scenario was 
also introduced for SN 1993 J by Fransson et al. (1996) to 
resolve a similar controversy between the theoretical high 
wind ionisation and a lower ionisation implied by the CS 
spectral lines. 

In Fassia et al. (2001) we suggested that the fast CS gas 
in SN 1998S resulted from the progenitor going through an 
RSG phase and then a BSC phase shortly before it exploded. 
'Mid-CSM' and 'Outer-CSM' zones were defined correspond- 
ing, respectively, to a tenuous, fast CSM wind from the BSC, 
and a dense, slow CSM wind from the RSG wind. The BSC 
wind would have 'piled up' against the inner boundary of 
the RSG wind, producing a negative velocity gradient. The 
SN shock would first engulf the faster more tenuous part of 
the BSG wind thus accounting for the observed decrease in 
the absorption blueshift with time. The abrupt appearance 
of the strong broad lines by about day 45 would be due to 
the impact of the ejecta with the denser RSG wind. Thus, 
qualitatively at least, this model can account for much of the 
observed spectral behaviour. However, an important result 
of the spectral line modelling described in Section 2 is that, 
in terms of the photospheric radius, the radius, R2, of the 
interface between the accelerated (fast) and unaccelerated 
(slow) wind remains about constant between days 23 and 
42 (Table 2). Yet, in the RSG/BSG wind explanation, R2 
decreases from 2.1 on day 23 to 1.3 on day 42. But Model 4 
(Table 1 and Fig. 2) shows that changing R2 from 1.3 to 
just 1.5 produces a much poorer reproduction of the ob- 
served profiles. This appears to present a difficulty for the 
BSG/RSG wind scenario. Moreover, the acceleration model 
can also account for the sudden appearance of the strong, 
broad lines simply as a result of the CDS becoming optically 
thin around this epoch (Chugai 2001). In other words, there 
may not be a need for a division of the wind into 'Mid-CSM' 
and 'Outer-CSM' components. All this leads us to favour 
the proposition that the fast CS gas is not a pre-existing 
flow, but instead is created by the ongoing acceleration of 
the slow wind in the immediate vicinity of the photosphere. 
The question of whether or not the BSG/RSG explanation 
can survive quantitative modelling is beyond the scope of 
this paper. 

It is instructive to review briefly the other known cases 
of type Iln and type IIL supernovae which have high ve- 
locity CS spectral features. The first observational evidence 
for a phenomenon of this kind was provided by the lUE 
spectra of the type IIL SN 1979C. This revealed UV lines 
of C IV and Si IV indicating a gas velocity in the range 
1000 - 4000 km s"^ (Panagia et al. 1980). Bearing in mind 
the photometric and spectral similarities between SN 1998S 
and SN 1979C, the fast (CS) material in SN 1979C may also 
be associated with shocked wind clumps. Originally, the UV 
lines in SN 1979C were attributed to pre-shock radiative ac- 
celeration of the slow wind (Chevalier 1981), or to emission 
from X-ray heated ejecta (Fransson 1984). 

Another case is the type lib SN 1993J where Ha 
P Cygni profiles suggested a gas velocity of ~ 1000 km s~^ 
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at 2-3 days post-explosion (Gumming et al. 1994). A simi- 
lar velocity was indicated by the N V 1240 A emission line 
(Sonneborn et al. 1994). Fransson et al. (1996) suggested ra- 
diative acceleration of the slow wind as a mechanism for the 
origin of this fast CS gas, but noted that Thomson scatter- 
ing could be insufficient. It would be interesting to determine 
whether or not line absorption radiation pressure could pro- 
vide the required enhancement of the radiative force in this 
case. Or do we need alternative mechanisms such as a cosmic 
ray precursor or shocked wind clouds? 

In the type Iln SN 1984E, spectral features are at- 
tributed to a superwind. Velocities of < 1000 km s~^ are 
indicated by the emission components of Balmer lines, while 
the blue edge of the Hq absorption indicates a velocity as 
high as 3000 km s~^ (Dopita et al. 1984), i.e. substantially 
higher than in the case of SN 1998S. The low signal-to-noise 
ratio of the SN 1984E data leaves some doubts about the 
precise interpretation of the line profile; in any event a ve- 
locity of 3000 km s~^ at the observed epochs is probably 
beyond the capability of radiative acceleration. 

The type Iln SN 1994W exhibited P Cygni lines which 
indicated CS gas velocities of 10^ km s~^. At first sight, 
therefore, we might take this SN to be similar to SN 1998S. 
However, the amazing persistence of the absorption line 
velocities for about two months (SoUerman, Gumming & 
Lundqvist 1998) is quite unlike the velocity evolution ob- 
served in the fast CS gas features of SN 1998S. A pre-existing 
CS shell (wind) with a constant velocity of ~ 1000 km s~^ 
is the more likely explanation for the fast CS gas around 
SN 1994W (SoUerman et al. 1998). 

It may well be that in different SNe Iln we see two kinds 
of fast CS matter. One type is the result of acceleration of 
the slow progenitor wind at the ejecta-wind interface (e.g. 
SN 1998S), while the other type is produced by the progen- 
itor before it explodes (e.g. SN 1994W). The origin of such 
different modes of heavy mass loss in the final stages of a 
star's life remains quite uncertain, with the conceivable in- 
volvement of stars from both the low mass (~ 10 Mq) and 
high mass (> 20 M©) ends of the massive star range. By 
studying a sufficiently large sample of SNe Iln using high- 
resolution spectroscopy together with other observations in 
the X-ray, optical and radio bands, we shall acquire clues 
about the complicated physics of the ejecta-wind interac- 
tion, the origin of the dense CS matter, and the pre-SN Iln 
genesis. 
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